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Abstract. We present the first determination of the proper motion of the neutron star low mass X-ray binary Cen X-4 measured 
from relative astrometry of the secondary star using optical images at different epochs. We determine the Galactic space velocity 
components of the system and find them to be significantly different from the mean values that characterize the kinematics of 
stars belonging to the halo, and the thin and the thick disc of the Galaxy. The high metallicity of the secondary star of the 
system rules out a halo origin and indicates that the system probably originated in the Galactic disc. A statistical analysis of 
the galactocentric motion revealed that this binary moves in a highly eccentric (e m 0.85 ± 0.1) orbit with an inclination of 
=* 110° to the Galactic plane. The large Galactic space velocity components strongly support that a high natal kick as a result 
of a supernova explosion could have propelled the system into such an orbit from a birth place in the Galactic disc. The high 
Li abundance in the secondary, comparable to that of stars in star forming regions and young stellar clusters like the Pleiades, 
may suggest a relatively recent formation of the system. Following the orbit backwards in time, we found that the system could 
have been in the inner regions of the Galactic disc ~ 100-200 Myr ago. The neutron star might have formed at that moment. 
However, we cannot rule out the possibility that the system formed at a much earlier time if a Li production mechanism exists 
in this LMXB. 
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1. Introduction 

Cen X-4 (V822 Cen) was discovered by the Vela 5B satellite 
during an X-ray outburst in 1969 (Conner et al. 1969). During 
the decay of a second outburst in 1979, a type I X-ray burst 
was observed, which indicated that the X-ray source is a neu- 
tron star (Matsuoka et al. 1980). This system is a prototype low 
mass X-ray binary (LMXB) consisting of a compact object of 
M NS = 0.5-2. 1 M Q (Shahbaz et al. 1993) and a donor star with 
a mass in the range 0.04 M Q < M 2 < 0.58 M Q (Torres et al. 
2002). 

The high Galactic latitude (b ~ 24°) of the system coupled 
with its distance (~ 1.2 + 0.3 kpc, Chevalier et al. 1989) places 
the object at the Galactic height of ~0.4-0.6 kpc and might 
suggest that it belongs to the Galactic halo. However, its rela- 
tively large centre-of-mass radial velocity (y = 191 + 1 kms~', 
Casares et al. 2005, in preparation) opens the possibility that 
the system received a natal kick as a result of the supernova 
explosion that generated the compact object and it might come 
from the Galactic disc. Such kicks have been proposed to ex- 
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plain the large transverse motions of neutron stars on the plane 
of the sky (e.g. Lyne & Lorimer 1994). 

The kinematic properties of other similar systems at high 
Galactic latitude have been studied by Mirabel et al. (2001, 
XTEJ1 118+40) and Mirabel & Rodrigues (2003, ScoX-1) 
by measuring their proper motions. These authors argued that 
XTE Jl 118+40 and ScoX-1 have kinematics consistent with 
the most ancient stars and globular clusters in the inner Galactic 
halo. However, according to the the galactocentric orbits of 
these objects, an origin in the bulge or disc of the Galaxy can- 
not be ruled out. In fact, Gualandris et al. (2004) have recently 
suggested that the black hole binary XTE J 1 1 1 8 +40 could have 
formed in the Galactic disc ~2-5 Gyr ago. 

The kinematic properties of Cen X-4 can help us to under- 
stand its origin. Here, we report the first determination of the 
proper motion of this system via relative astrometry of the sec- 
ondary star using optical images. This, together with the radial 
velocity and the distance, allows us to determine the kinematic 
parameters and study the galactocentric orbit of this important 
LMXB. 
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2. Proper motion 

Accurate proper motion measurements of LMXBs have been 
obtained using very long baseline array (VLBA) or very 
long baseline interferometry (VLB1) techniques at radio wave- 
lengths and/or from optical observations of these systems 
at HST or large telescopes. Unfortunately, there are no re- 
ported observations of Cen X-4 with either of these facilities. 
However, several campaigns have carried out observations of 
this target using optical ground-based telescopes. We selected 
two images among those with better pixel scales: a plate from 
the Second Epoch Survey 1 of the Southern Sky observed in 
May 1994 in the Anglo- Australian Observatory (AAO) with 
the UK Schmidt Telescope with a pixel scale of 1 arcsec, and 
a second image with a pixel scale of 0.238 arcsec obtained 
in May 2002 with the Wide Field Imager at the European 
Southern Observatory (ESO), Observatorio La Silla, using the 
Cassegrain focus of the 2.2 m MPG/ESO telescope. 
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Fig. 1. Distribution of shifts in right ascension and declination 
for the 39 reference stars and the target (wide cross), Cen X-4. 
The error bars give the total error (see §2). 

The proper motion was determined using several tasks 
within IRAF. First, we measured the pixel coordinates of 40 
stars including the target, Cen X-4, in both images (with the 
task imcentroid). Second, the 1994 image was registered on 
the 2002 one, by applying the coordinate transformation (us- 
ing the tasks geomap and geoxytrans) between both frames 
of reference objects. We selected as reference objects those 
stars whose transformation residuals were lower than 2cr above 
the rms of the transformation (first criterion). Third, we mea- 
sured the relative shifts of all the stars between their positions 
corresponding to the 1994 and 2002 images in the reference 
frame of the 2002 image. We then rejected those stars whose 
proper motions were higher than 3<x above the standard devia- 
tion of the motions of all the stars (second criterion). The whole 
procedure was iterated until we could match the reference ob- 

1 Plates from this survey, also called Second Palomar Sky Survey 
(POSSII), have been digitized and compressed by the ST Scl. The dig- 
itized images are copyright (c) 1993-2000 by the Anglo- Australian 
Observatory Board and are distributed herein by agreement. All Rights 
Reserved. Produced under Contract No. NAS5-2555 with the National 
Aeronautics and Space Administration. 



jects with both criteria unambiguously. We aligned (using the 
tasks ccmap and cctran) the 2002 frame in right ascension and 
declination with these reference objects using the coordinates 
given by 2MASS 2 . In order to correct for geometrical distor- 
sions with geomap and ccmap, we selected general geometry 
transformations up to order 2. These take into account rota- 
tion, pixel/celestial coordinate scale and shift between the ref- 
erence frames of both images. Finally, we computed the rela- 
tive displacements of our target and estimated the error from 
the average relative motions of the reference stars (see Fig. 
1). The computed proper motion is fi a = 11 + 10 mas/yr and 
Us = -56 +10 mas/yr, which leads to an overall yearly motion 
of pL = -57 ±12 mas/yr along a position angle with respect to 
the North Celestial Pole of 100° ± 12°. At a distance d (in kpc), 
this proper motion corresponds to a transverse velocity on the 
plane of the sky of (270 ± 56)d km s _1 . 

3. Distance to Cen X-4 

In order to determine the space velocity of the system from the 
proper motion measurement we need to know its distance to 
the Sun. The optical properties of the secondary star can pro- 
vide an estimate of the distance, although we found several 
discrepancies among the values obtained using different photo- 
metric filters. The mass range estimated for the secondary star 
of 0.04 < M 2 /M e < 0.58 (Torres et al. 2002) leads to a stellar 
radius of 0.5 < R2/R0 < 1-2 by assuming that the secondary 
fills its Roche lobe. This radius, together with our spectroscopic 
estimate of the effective temperature, T e s = 4500 ± 100 K 
(Gonzalez Hernandez et al. 2005), provides an intrinsic bolo- 
metric luminosity of 0.095 < L2/L© < 0.371. The apparent 
magnitude is my ~ 18.2 (van Paradijs et al. 1980) and the 
colour excess, E B - V = 0.1 ± 0.05 (Blair et al. 1984). Using 
the bolometric corrections (Bessel et al. 1998) and the effec- 
tive temperature, surface gravity and veiling from the accre- 
tion disc derived by Gonzalez Hernandez et al. (2005, submit- 
ted), we found 1.5 < d mv (kpc) < 3.6. Note that Cen X-4 is 
highly variable in quiescence, with reported magnitudes my ~ 
18-18.7 (Chevalier et al. 1989; Cowley et al. 1988). From the 
m R ~ 17.5 and m H ~ 15.1 (Shahbaz et al. 1993) and assum- 
ing A R /A V ~ 0.83 and A H /A V ~ 0.21 (Schaifers et al. 1982), 
we obtained 1.2 < d mR (kpc) < 2.9 and 0.9 < d mH (kpc) < 2.3. 
The distance obtained from the V and R filters in particular are 
rather uncertain because our veiling estimates and the assumed 
apparent magnitudes were determined at different epochs and 
are also affected by the contribution of a contaminant star (with 
m v < 20.5) to the total flux (Casares et al. 2005, in prepara- 
tion). 

A more confident value of the distance can be calculated 
from the peak fluxes of the X-ray precursor observed in the 
1969 outburst (Belian et al. 1972) and the X-ray burst detected 
during the 1979 outburst (Matsuoka et al. 1980). From the ob- 
served fluxes of ~ 1.4 x 10~ 6 erg crrT 2 s _1 and ~ 1 x 10~ 6 

2 The Two Micron All Sky Survey is a joint project of the 
University of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science 
Foundation. 
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Fig. 2. Distribution of perigalacticon and apogalacticon distances (upper left panel), maximum vertical distance (upper right 
panel), eccentricities (lower left panel) and orbital inclinations (lower right panel). Distributions are normalized to the total 
number of orbits in our calculus (N t = 5000). 



erg crrT 2 s _1 (for 1969 and 1979, respectively) and by assum- 
ing that they correspond to X-ray luminosities at the Eddington 
limit for a canonical neutron star (Z-Edd ~ 1.6-2.7 xlO 38 erg 
s" 1 ), Chevalier et al. (1989) derived a distance of 1.2 + 0.3 kpc. 
This measurement can also be estimated by considering recent 
average peak luminosities of photospheric radius expansion 
bursts observed in globular clusters, which is (2-3.8) xlO 38 erg 
s _1 (Jonker & Nelemans 2004, and references therein). From 
these peak luminosities and the observed fluxes we obtained a 
distance of d - 1.4 + 0.3 kpc. 

The space velocity calculated according to the derived 
proper motion exceeds the Galactic escape velocity (around 
500 km s _1 in the solar neighbourhood) if the system is fur- 
ther than 1.7 kpc from the Sun (see §4). Therefore, we assume 
hereafter a distance to Cen X-4 of 0.9 < J(kpc) < 1 .7 which is 
where both the X-ray and optical distance measurements over- 
lap. 

4. Space velocity 

From the values of position, distance, proper motion and ra- 
dial velocity, the velocity components J/ L sr, Vlsr an d Wlsr 
relative to the Local Standard of Rest can be calculated us- 
ing, for example, Johnson & Soderblom's (1987) equations 
of transformation. Assuming that the Sun moves with compo- 
nents (£/ 0> lsr, V 0> lsr, W g ,lsr) = (10., 5.2, 7.2) km s -1 (Dehnen 
& Binney 1998) relative to the Local Standard of Rest, from 



the two extreme values of the distance we obtained: for d = 0.9 
kpc, (t/ LSR , Vlsr, W lsr ) = (204 ± 21,-196 ± 35,-122 ± 38) 
km s-\ and for d = 1.7 kpc, (U hSR , V LS r, W lsr ) = (239 ± 
40, -303 ± 67, -305 ± 71) km s _1 . This implies that the source 
moves towards the Galactic Centre and has a 3cr vertical mo- 
tion towards the Galactic plane. These two sets of values are 
rather different from the mean values of U, V and W that carac- 
terize the kinematics of stars belonging to the halo, and the thin 
and thick discs of the Galaxy (Chiba & Beers 2000; Soubiran 
et al. 2003). The system velocities would be v sys = 308 + 51 
kms" 1 and v sys = 493 + 105 kms" 1 for d = 0.9 kpc and d = 1.7 
kpc respectively. 

5. Galactocentric orbit 

We integrated the orbit of Cen X-4 based on the heliocentric 
position and kinematic data listed in Table 3. In order to con- 
vert heliocentric positions and velocities into coordinates in the 
Galaxy frame we placed the Sun at (X Q , Y G ,Z ) = (-8,0,0) 
kpc, with a velocity r Q = (U Q ,V Q ,W Q ) = (10.0,225.2,7.2) 
(Binney & Merrifield 1998; Dehnen & Binney 1998) and used 
the equation 

r = r Q + v r (cos b cos /, cos b sin /, sin b) 
+r s /diK(- sin cos 0) 
+r s ^bK(- sin b cos /, - sin b sin /, cos b) 
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Fig. 3. Galactocentric orbits of Cen X-4 viewed from above the Galactic plane (left panels) and edge on (centre and right panels). 
The upper panels show the orbit for a heliocentric distance of d = 0.9 kpc and the lower panels for d = 1.3 kpc. In the left 
panels, the solid line represents the section of the orbit in the Southern Galactic Hemisphere, whereas the dashed line shows the 
section of the orbit in the Northern Galactic Hemisphere. Both orbits were calculated by integrating backwards over 460 Myr. 
The circular orbit corresponds to the Sun. Open and filled circles show the current positions of the Sun and Cen X-4 respectively. 



where r s is the heliocentric distance, k 4.74, a conversion 
factor from (kpc mas/yr) to (km s _1 ), v r , the heliocentric ra- 
dial velocity and Hi,fib> proper motions in the longitudinal and 
latitudinal direction, respectively. 

We assumed a three-component Galaxy model that consists 
of: (i) a Miyamoto & Nagai (1975) disc, (ii) a Hernquist bulge 
(Hernquist 1990) and (iii) a spherical logarithmic halo. The 
Galaxy parameters can be found, for example, in Law et al. 
(2004). 

Unfortunately, the heliocentric distance and velocity of 
Cen X-4 are poorly determined. Values with these large errors 
represent a large set of possible orbits that must be statistically 
treated. We calculated an initial set of 5000 orbits using Monte 
Carlo simulations, where the heliocentric distance is randomly 
chosen within the range given above. Proper motions were dis- 
tributed following a Gaussian with variance equal to the mea- 
sured error of 10 mas/yr and mean value (///) = -22 mas/yr 
and (fib) = -53 mas/yr. Subsequently, we integrated each orbit 
1 Gyr backwards in time. In Fig. 2 we show the orbital parame- 
ters that result from the distribution of initial positions and ve- 
locities. The orbit of Cen X-4 appears totally inconsistent with 
typical orbits of Milky Way stars. This system has a perigalac- 
ticon of r p ~ 1 - 4 kpc, whereas the apogalacticon occurs at 
r a ~ 14-16 kpc from the Galactic Centre. Although r p and r a 
appear to be fairly sensitive to the initial position and velocity, 
the resulting eccentricity, e = (r a - r p )/(r a + r p ) = 0.85 + 0.10, 
barely depends on the errors from heliocentric distance and 



proper motion measurements. This eccentricity is considerably 
higher than that of disc stars (e 0.5-0.1). The orbital inclina- 
tion from the disc plane also shows a well defined distribution 
peak at ; 110° and determines its orbital sense of motion to 
be retrograde for most of the distances considered in the dis- 
tribution. This orbital inclination is considerably higher than 
for Milky Way stars (i.e. i Q = 4°) and leads to high vertical 
motions. As we can see in the upper right panel, Cen X-4 will 
move (or may have moved) in halo regions, with a most likely 
maximum vertical distance as large as 12-16 kpc. 

In Fig. 3 we have represented the galactocentric orbits of 
Cen X-4 for two distances, 0.9 and 1.3 kpc. These orbits were 
calculated by integrating backwards in time for 460 Myr and 
both orbits are highly eccentric. The orbit for d = 0.9 kpc 
shows a prograde sense of motion, whereas for d = 1.3 kpc 
the system describes a retrograde orbit. These orbits intersect 
the disc plane for the first time (integrating backwards in time) 
148 and 243 Myr ago for d = 0.9 kpc and d = 1.3 kpc, respec- 
tively. 

6. Pinning down the formation region of Cen X-4 

Two main scenarios have been proposed to explain the current 
kinematic properties of LMXBs observed at high Galactic lati- 
tudes: (i) these binary systems originated in the Galactic bulge 
or in the disc and received a natal kick as a result of black 
hole/neutron star formation, or (ii) they formed in the Galactic 
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Fig. 4. Upper panel: First crossing of the axisymmetry plane when integrating the orbit of Cen X-4 backwards in time. The 
circle indicates its present position, and its diameter corresponds to the uncertainty in the distance estimate. Lower left panel: 
Distribution of times when the orbit (integrated back in time) crosses the disc plane for the first time. Lower right panel: distri- 
bution of energies per unit of mass gained by Cen X-4 in comparison with circular orbits at the first-crossing position. The solid 
line represents the definition for the energy E given in the text whereas the dotted and dashed lines show E' for M He = 4 M Q and 
M 2 = 0.8 M Q for M He = 3 M and M 2 = 0.5 M respectively. 



halo either in dense globular clusters or as isolated ancient pri- 
mordial binaries. Mirabel & Rodriguez (2003) argued that the 
LMXB Sco X-l most probably was formed in the core of a 
globular cluster in a close encounter of the compact object with 
a single star or with a binary (see Verbunt 2003 for a review). 

Halo stars and globular cluster stars are generally metal- 
poor (with heavy elements between 10 and 1000 times less 
abundant than in the Sun). If the Cen X-4 system had origi- 
nated in the halo, the secondary star would have a low metal- 
licity. Recently, the iron content 3 of the secondary star has 
been determined at roughly twice solar ([Fe/H] = 0.23 + 0.1, 
Gonzalez Hernandez et al. 2005). This high metal content, to- 
gether with measured element abundance ratios in the compan- 
ion star, makes it very unlikely that the system was formed 
in the Galactic halo, either in a globular cluster or as an an- 
cient primordial binary. Most probably, the system formed in 
the Galactic disc and was expelled from its birth place as a 
consequence of a supernova (SN hereafter) explosion of the 
massive progenitor of the neutron star. 

An origin in the Galactic disc requires that the system ob- 
tained enough energy to be expelled, most probably as a conse- 
quence of a natal kick in a SN explosion. In the lower right 

3 [Fe/H] = log[iV(Fe)/Ar(H)L - log[JV(Fe)/JV(H)] . 



panel of Fig. 4 we show an estimate of the energy per unit 
of mass gained by the binary system due to a supernova at a 
first crossing. We define this energy as E = E orb - E cilc , where 
E C i K is the energy per unit of mass of a star following a cir- 
cular orbit in the region where Cen X-4 crosses the Galactic 
plane. The resulting distribution is close to a Gaussian with 
(E) - (8 + 1) x 10 4 (kms _1 ) 2 . This energy must be lower than the 
difference between the total energy of the system after and be- 
fore the SN explosion if the binary components remain bound. 
Brandt & Podsiadlowski (1995) investigated the effects of high 
SN kick velocities on the orbital parameters of post-SN neutron 
star binaries. Using their formalism we compared the pre-SN, 
E vk -sn, and post-SN, £p OSt -sN, total energies per unit of mass 
of the system after undergoing SN explosions. We defined this 
energy balance as E' = £po S t-SN(M T> p OSt _ SN /M T ,p re -SN)-£pre-SN 
where M TiPOSt _sN and M TiPre _sN are the total masses of the 
system after and before the SN event respectively. We found 
that energy released in a SN of He core progenitor of mass 
M He = 4 M G and a secondary star of M 2 = 0.8 M is 
E' ~ 1 - 3 x 10 4 (km s" 1 ) 2 , whereas for M He = 3 M and 
M 2 = 0.5 M , E' ~ 1 - 2 x 10 4 (km s" 1 ) 2 . These energies are 
consistent with the energy required to launch Cen X-4 into an 
orbit like this from a birth place in the Galactic disc (see Fig. 
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4). A detailed study of the distribution of the kick directions 
with respect to the orbital motions of each component shows 
that it strongly depends on the masses of the different compo- 
nents of the system before and after the SN explosion. Thus 
the kick direction cannot randomly be chosen and is more re- 
stricted when the differences between the masses of both com- 
ponents are higher. 

The abundance of lithium in the secondary star of Cen X-4 
(log e(Li)crE = 3.06+0.29) is close to the cosmic Li abundance 
in the Galactic disc (Martin et al. 1994a; Gonzalez Hernandez 
et al. 2005) and thus substantially higher than in stars of com- 
parable mass of the young Pleiades cluster whose age is ~ 120 
Myr (Martin et al. 1998, Stauffer, Schultz & Kirkpatrick 1998). 
This might indicate that the system is relatively young since 
Li is severely depleted in the atmospheres of such stars during 
pre-main sequence and main sequence evolution. In order to 
explore the possibility of a recent formation of Cen X-4 in the 
Galactic disc, we calculated the positions where the orbit inter- 
sects the Galactic plane for the first time (integrating backwards 
in time). The result is plotted in the upper panel of Fig. 4, which 
shows that the orbit crosses the inner region of the Milky Way 
in the range 0—4 kpc from the Galactic Centre. The distribu- 
tion of times (lower left panel) indicates that the system was 
probably in the Galactic disc ~ 100-200 Myr ago. The high Li 
abundance of the secondary star might be explained if the sys- 
tem had formed at that moment. However, it has been proposed 
that a Li production mechanism via spallation reactions could 
be responsible for the high Li abundances of LMXBs (Martin 
et al. 1994b; Yi & Narayan 1997; Guessoum & Kazanas 1999). 
No direct evidence for such a mechanism has been found yet in 
Cen X-4 or in any other LMXB. Nevertheless, we cannot rule 
out that such a mechanism has been active in Cen X-4 and, 
thus, that the system could be much older. 

7. Conclusions 

We have used optical relative astrometry of the secondary star 
of the neutron star binary Cen X-4 to measure its proper mo- 
tion. This proper motion, together with the estimated distance 
and radial velocity and its celestial coordinates, allows us to 
determine the space velocity components. These appear to be 
rather different from the mean values of U, V and W for stars 
belonging to the halo, the thin and the thick disc of the Galaxy. 
The system is falling into the Galactic plane with a large verti- 
cal velocity > 100 km s _1 . 

The high metallicity ([Fe/H] = 0.23 + 0.1, Gonzalez 
Hernandez et al. 2005) and the element abundance ratios of 
the companion star in Cen X-4 makes it very unlikely that the 
system formed in the Galactic halo. Most probably, the system 
originated in the Galactic disc and was propelled into its cur- 
rent position by a natal kick in a supernova explosion. 

We have explored the galactocentric orbit of Cen X-4 
through a statistical analysis that takes into account the un- 
certainty in the proper motion and the distance. The orbital 
inclination from the Galactic plane is 110° and determines 
its orbital sense of motion to be retrograde in a highly eccen- 
tric (e =a 0.85 + 0.1) orbit. Following this orbit backwards, we 
have found that the system was probably at the Galactic disc 



~ 100-200 Myr ago. The formation of the neutron star at that 
time would imply a relatively young system and may provide 
a natural explanation for the observed high Li abundance of 
the secondary star. However, if there is an active Li production 
mechanism in Cen X-4, the Li abundance cannot be directly in- 
terpreted as an indication of youth, and it cannot be discarded 
that the system formed at much earlier times. 
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